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20>  ABSTRACT  (Cmatkmtm  mt»  rartrw  iMi  tt  «**d  Identity  by  block  ntmtbmr) 

~  Empirical  equations  for  the  drift  velocities  of  electrons  and  holes  in  silicon  are  given  as  a  func¬ 
tion  of  electric  field,  temperature,  and  doping  density.  A  single  equation,  which  Is  valid  above 
room  temperature,  results  from  the  inverse  dependence  of  the  saturation  velocity  upon  the 
square  root  of  the  temperature. 
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1.  Introduction 


(3) 


Simple  empirical  equations  for  the  drift  veloc¬ 
ities  of  electrons  and  holes  as  a  function  of  electric 
fields,  temperature,  and  doping  levels  are  needed 
for  device  simulation.  A  review  of  charge  transport 
properties  of  silicon,  including  phenomenological 
expressions,  was  given  by  Jacoboni  et  alT  How¬ 
ever,  none  of  their  expressions  included  all  three  of 
the  variables  of  interest.  Schwarz  and  Russek2 
presented  semi-empirical  equations  for  electrons 
in  silicon,  but  their  expressions  are  not  simple  to 
use.  The  well-known  empirical  equation  of  Schar- 
fetter  and  Gummel3  is  for  room-temperature  only. 

2.  Derivation  of  Equation 

The  equation  of  Scharfetter  and  Gummel3  for 
the  drift  velocity,  vd,  as  a  function  of  field,  E,  and 
doping  level,  N,  is 

vd  =  MoE/({l  +N/[(N/S)  +  NrJ 

+  (E/A)2/[(E/A)  +  F]  +  (E/B)2»  ,  (1 ) 

where  h0  is  the  low-field  mobility  and  all  other  sym¬ 
bols  are  fitting  parameters.  The  temperature  de¬ 
pendence  of  nD  is  generally  given  as 

f'o-CT-r,  (2) 

where  T  is  temperature  and  C  and  y  are  constants. 
The  power  dependence  varies  according  to  the 
investigator,  but  is  usually  around  2.5  for  both  holes 
and  electrons. 

The  temperature  variation  of  the  saturation  ve¬ 
locity,  vs,  for  electrons  was  measured  by  Duh  and 
Moll. 4  Their  data  are  plotted  on  a  log-log  scale  in 
figure  1 .  It  is  seen  that  there  is  a  close  fit  to  a  power 
law,  i.e., 
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vs«k/p\ 

in  the  range  above  room  temperature. 

An  empirical  expression  for  vs  is  given  by 
Jacoboni  et  alf  as 

vs  =  v*/[l  +  G  exp  (TIQ)] ,  (4) 

where  v*  *  2.4  *  107  cm/s,  G  =  0.8,  and  0  =  600  K. 
This  expression  is  also  plotted  in  figure  1 . 

The  saturation  velocity  from  (1)  is  found  to  be 

Vs  *  Mob  ■  (5) 

By  equating  (3)  and  (5)  we  obtain,  from  (2), 

B  =  DTr~Vi ,  (6) 

where  D  *  K/C.  Finally,  using  (2)  and  (6)  in  (1)  we 
obtain  our  desired  expression, 

vd  *  (CT  TE)/[  1  +  N/[(N/S)Nr] 

+  [(E/A)2]/[(E/A)  +  F]  +  (E/DT_Va)2],/l  (7) 


Figure  1.  Saturation  velocity  as  function  of  T.  Points  are 
measured  data  of  Duh  and  Moll.4  Dashed  line  Is  from  equa¬ 
tion  (4)  and  solid  line  from  equation  (3). 


3.  Results 

Plots  of  vd  from  equation  (7)  for  electrons  in  sili¬ 
con  are  given  in  figure  2  for  N  *  0  at  300,  400,  500, 


DRIFT  VELOCITY  (cm/t) 


and  600  K.  The  parameters  used  are  listed  in  ta¬ 
ble  1.  Also  plotted  in  figure  2  are  600-K  curves  for 
doping  levels  of  1  *  1017,  1  x  1018,  and  1  x  1019 
cm-3.  All  curves  are  seen  to  be  smooth.  Room- 
temperature  plots  of  vd  as  a  function  of  E  for  vari¬ 
ous  doping  levels  may  be  found  in  Jacoboni  et  al.? 
Plots  of  the  temperature  dependence  of  vd  at  fixed 
fields  show  a  gradual  decrease  in  the  exponent  of  T 
from  y  at  low  fields  to  1/2  at  high  fields.  Plots  have 
been  made  of  hole  drift  velocities  from  equation  (7) 
and  from  the  parameters  of  table  1 ,  and  they  also 
show  smooth  variations. 


Figure  2.  Electron  drift  velocity  from  equation  (7)  with 
parameters  of  table  1.  Solid  curves  are  for  N  «  0  and  con¬ 
stant  temperatures  (K)  as  labeled.  Dashed  curves  are  for 
T  =  600  K  and  doping  densities  (cm-3)  as  labeled. 


Table  1.  Parameters  used  in  equation  (7)  for 
calculating  drift  velocities 


Parameter 

Electrons 

Holes 

Units 

A 

3  5  *  1 03 

6.1  x  io3 

V/cm 

c 

1  43  *  109 

1  35  x io8 

cm2K7V~'s-1 

D 

1  29  * 10_1 

1.37 

cmK  l/,-7s-1 

F 

88 

16 

— 

Nr 

3  x  1016 

4  x  1016 

cm-3 

S 

350 

81 

— 

Y 

2  42 

2  2 

— 

'C  Jacoboni,  C  Canah,  G  Ottaviam,  and  A  Alberigi  Quaranta. 
A  Review  of  Some  Charge  Transport  Properties  of  Silicon.  Solid- 
State  Electron .  20(1977).  77  89 


Cook  and  Frey5  have  measured  the  high-field 
electron  drift  velocity  in  silicon-on-sapphire  (SOS) 
films  doped  to  N  =  2  x  1017  cm-3.  Their  data  are 
plotted  in  figure  3.  Also  plotted  in  this  figure  are  cal¬ 
culations  based  on  equation  (7)  and  the  parameters 
of  table  1 .  The  low-field  mobility  was  reduced  by  50 
percent  in  a  second  calculation,  also  plotted  in 
figure  3.  To  show  the  flexibility  of  (7),  the  parameter 
F  was  set  at  <*>  (i.e.,  the  third  term  of  the  denomina¬ 
tor  was  set  to  zero),  and  the  lower  mobility  calcula¬ 
tion  was  repeated.  The  fit  to  the  experimental  data 
is  quite  good.  Schwarz  and  Russek2  predict  the 
steeper  rise  in  the  velocity-field  profile  in  doped 
silicon.  When  further  data  are  available,  the  param¬ 
eters  A  and/or  F  could  perhaps  be  made  dependent 
upon  the  doping  level. 


Figure  3.  Electron  drift  velocity  as  function  of  field  for 
N»2*10-17  cm-3  and  room  temperature.  Experimental 
points  with  error  bars  are  from  Cook  and  Frey.5  Upper 
dashed  curve  Is  from  equation  (7)  with  table  1  parameters. 
Lower  dashed  curve  is  calculated  with  lower  mobility,  but 
same  saturation  velocity.  Solid  curve  shows  latter  mobility 
with  third  term  of  denominator  of  equation  (7)  omitted. 
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